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ABSTRACT 

The standard solar model was so reliable that it could predict the existence 
of the massive neutrino. Helioseismology measurements were so precise that 
they could determine the depth of the convection zone. This agreement between 
theory and observation was the envy of all astrophysics - until recently when so- 
phisticated three-dimensional hydrodynamic calculations of the solar atmosphere 
reduced the metal content by a factor of almost two. Antia & Basu (2005) sug- 
gested that a higher value of the solar neon abundance, A Ne /A Q = 0.52, would 
resolve this controversy. Drake & Testa (2005) presented strong evidence in favor 
of this idea from a sample of 21 Chandra stars with enhanced values of the neon 
abundance, A Ne /A Q = 0.41. In this paper, we have analyzed solar active region 
spectra from the archive of the Flat Crystal Spectrometer on Solar Maximum 
Mission, a NASA mission from the 1980s, as well as full-Sun spectra from the 
pioneering days of X-ray astronomy in the 1960s. These data seem consistent 
with the standard neon-to-oxygen abundance value, A Ne /A Q = 0.15 (Grevesse 
& Sauval 1998). If these results prove to be correct, than the enhanced-neon 
hypothesis will not resolve the current controversy. 

Subject headings: stars: abundances — stars: coronae — Sun: abundances - 
Sun: corona — Sun: X-rays, gamma rays 

1. Introduction 

Scientific progress is made as theory evolves to predict observables and observations 
confirm theoretical predictions. There was no better astronomical example of this progress 
than the agreement between the complex models and helioseismology measurements of the 
solar interior - that is, until recently. New 3-D hydrodynamic models (Allende Prieto, Lam- 
bert & Asplund 2001; Asplund et al. 2000) which account for temperature inhomogeneities, 
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convection and the associated velocity fields, and relax the assumption of local thermody- 
namic equilibrium have been used to compute atomic level populations. When these results 
are applied to solar absorption lines, the photospheric abundances of light elements such as 
C, N, O, and Ne all must be reduced by 25-35% (Asplund et al. 2004). But these elements 
provide a major source of opacity for the solar interior, which determines the internal so- 
lar structure and the depth of the convection zone. The changes are significant enough to 
shatter the harmonious agreement between theory and observation that had come to exem- 
plify scientific progress. The discrepancies between the two are now far in excess of either 
the observational uncertainties or the model predictions (Bahcall, Serenelli & Pinsonneault 
2004). 

Antia & Basu (2005) and Bacall, Basu & Serenelli (2005) have suggested that one way 
to restore the disrupted harmony is to increase the solar abundance of neon. Where the C, N, 
and O abundances are determined from photospheric absorption lines, the neon abundance 
must be inferred from solar energetic particles and the spectra of hot stars and nebulae; it 
is often scaled to the oxygen abundance: A^ e jAo = 0.15 (Grevesse & Sauval 1998). If this 
value were 3.44 times higher, or A^ e j Ao = 0.52, the discord between the solar models and 
the helioseismology observations would disappear. 

To this end, Drake & Testa (2005) determined the Ne/O abundance ratios for a sample 
of 21 stars located within 100 pc of the Sun which had been observed with the High En- 
ergy Transmission Grating on the Chandra X-ray Observatory. These results were possible 
because the X-ray spectra of cool stars (including the Sun) feature prominent emission lines 
of highly ionized ions of both oxygen and neon. These authors correctly point out that the 
abundances of neon and oxygen may be different in the photospheres and coronae of these 
stars - they are both elements with high First Ionization Potential (FIP) - but the ratio 
of the abundances should be identical. Their Chandra results suggest that A Ne /A Q = 0.41, 
significantly higher than the accepted solar value, and high enough (assuming other un- 
certainties) to resolve the controversies between the solar models and the helioseismology 
observations. In this paper, we have reexamined X-ray spectra similar to those of Drake & 
Testa (2005), but for the Sun itself. 



2. Analysis 

The spectra analyzed here were obtained with the Flat Crystal Spectrometer (FCS - 
Acton et al. 1980) on Solar Maximum Mission, which operated from 1980-89. The instru- 
ment had a 15 arcsec field-of-view and could scan the soft X-ray resonance lines of prominent 
ions in the range of 1.5 A to 20.0 A with a spectral resolution of 0.015 at 15 A. These lines 
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are sensitive to a broad range of plasma temperatures (1.5-50 MK) and can be used to deter- 
mine temperature, emission measure, density, and relative elemental abundances. The FCS 
observing sequence used to accumulate these spectra began with a low-resolution image of 
the active region. The instrument was then pointed at the brightest pixel of that image and 
began scanning at a rate of about 0.01 A per second from 13 A to 19 A. Shorter wavelengths 
were scanned simultaneously in different channels, but all the lines used in this Letter are 
in this limited range. Covering the full wavelength range took approximately 10 minutes, 
and four to five of these spectroscopic scans were done during the daylight portion of the 
90-minute spacecraft orbit. 

These FCS spectra have been used in the past to study elemental abundances, but differ- 
ent science goals were driving those analyses: Strong, Lemen & Linford (1991) searched for 
statistically significant abundance variations; Schmelz (1993) discovered abundance anoma- 
lies in flaring solar plasma; Fludra & Schmelz (1995) found a set of absolute coronal abun- 
dances; and Schmelz et al. (1996) looked for deviations from the standard FIP effect. In 
none of these cases were the data used to compute the normal solar coronal Ne/O abun- 
dance ratio because, at the time, that value was known. Now that doubt has been cast on 
the value of this ratio, we felt it was important to examine the FCS data archive to see if 
the data could confirm the higher Ne/O abundance ratios proposed by Antia & Basu (2005) 
and observed by Drake & Testa (2005). 

The FCS archive contains 75 of the relevant spectral scans described above. Those taken 
during flares and long duration events were excluded from this sample to eliminate physical 
conditions which changed with time. We also excluded those spectra where the hot Fe XIX 
A13.52 line (peak formation temperature Log T = 6.8 MK) was visible. The presence of this 
line would indicate dynamic, rather than quiescent, plasma conditions. Of the remaining 
spectra we selected 20 with the best signal-to-noise ratio in order to be as confident as 
possible that the fluxes we have measured are dominated by the lines of interest. We will 
be using three Ne IX lines, three O VIII lines, and a Fe XVIII-to-Fe XVII line ratio as a 
temperature diagnostic. The wavelengths, peak formation temperatures, and spectroscopic 
notation for these transitions are listed in Table 1. The FCS data for these lines are listed in 
Table 2, which includes the date and time of the observation, the solar active region number, 
the electron temperature determined from the flux ratio of Fe XVIII A14.21 to Fe XVII 
A16.78, as well as the fluxes and uncertainties in photon units for the relevant neon and 
oxygen lines. 

If we assume that the plasma in the FCS field-of-view is isothermal, then the observed 
flux, F, of a soft X-ray emission line is proportional to the elemental abundance, A, and the 
contribution function, G(T). Then the Ne/O abundance ratio can be computed as follows: 
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The contribution functions were taken from version 4.02 of the CHIANTI Atomic Physics 
database (Dere et al. 1997; Young et al. 2003), and we have used the ionization fractions 
of Mazzotta et al. (1997). The G(T) functions for the different neon and oxygen lines used 
in this analysis are plotted in the left column of Figure 1 and the Go{T)-to-G Ne{T) ratios 
for different pairs of the neon and oxygen lines are plotted in the right column. We can see 
from these plots that the ratios have a significant temperature dependence. Drake & Testa 
(2005) added a fraction of the Ne X A12.13 flux to the Ne IX A13.45 flux to remove much of 
this temperature dependence: G^e = G^e ix + 0.15GAr e x- Unfortunately, however, this 
Ne X A12.13 line was not part of the FCS spectral scan. If we assume that the plasma in 
the FCS field-of-view is isothermal, we can compute the temperature from the flux ratio of 
Fe XVIII A14.21 to Fe XVII A16.78; these results are listed in Table 2. However, it is much 
more likely that the observed plasma is multi-thermal, with a Differential Emission Measure 
(DEM) typical of a quiescent active region (see, e.g., Brickhouse, Raymond & Smith 1995; 
Brosius et al. 1996; Schmelz et al. 2001). If this is true, then the plasma temperatures listed 
in Table 2 are the mean values weighted by the DEM itself and the contribution functions 
of the two iron lines. 

Unfortunately, it is not possible to construct a full DEM for the observed active region 
plasma using the FCS lines alone. (This is possible for a flare, however; see, e.g., Schmelz 
1993). The problem is that the lowest temperature lines in the data set are the O VIII 
lines listed in Table 1, with a peak formation temperature of Log T = 6.5. Typical DEM 
curves of quiescent active region plasma have significant contributions at lower temperatures. 
So, although the high-temperature end of the DEM curve would be well-constrained by 
observations of the Fe XVIII A14.21 line flux and upper limits of the Fe XIX A13.52 line 
flux, the FCS lines, unfortunately, do not provide a useful low-temperature constraint. 

We circumvent this problem in the following way: although it is not possible to determine 
the exact Ne/O abundance ratio without assuming the shape of the DEM curve, it is possible 
to determine an upper limit for this value. So, instead of searching for the precise Ne/O 
abundance ratio, we seek to determine if the FCS data are consistent with a Ne/O abundance 
ratio as high as A^ e /Ao = 0.52, the value proposed by Antia & Basu (2005) which would 
eliminate the discord between the solar models and the helioseismology observations, or even 
as high as A Ne /A Q = 0.41, the value obtained by Drake & Testa (2005) using the Chandra 
observations for a sample of 21 nearby stars. 

To calculate this upper limit, we examine the plots in Figure 1. In the temperature 
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range of the typical quiescent active region DEM (6.0 < Log T < 6.6), the Go(T)-to-GNe(T) 
ratios in Figure 1 all show a local maximum at Log T w 6.3. Therefore, the upper limit of 
the Ne/O abundance ratio will result if we assume that the plasma temperature is Log T 
6.3. This assumption is made in the analysis shown in Figure 2 where we have used Eq. 1 to 
calculate A^ e /Ao with F^ e /Fo from the FCS observations for the 20 spectra listed in Table 
2, and Go(6.3)/GAr e (6.3) from the CHIANTI database. Each panel of Figure 2 also shows 
the weighted mean (solid line), the la uncertainty (dotted line), A Ne /A Q = 0.41 (dot-dash 
line) from Drake & Testa (2005), and A Ne /A Q = 0.52 (dash line) from Antia & Basu (2005). 
The panels of Figure 2 show the following ratios (a) Ne A13.45-to-0 A16.01; (b) Ne A13.55- 
to-0 A18.97; and (c) Ne A13.77-to-0 A15.18. The weighted mean and its uncertainty are 
printed in the upper left corner of each plot. (Note: we show only three Ne/O ratios here, 
but all ratio possibilities gave similar abundance results.) 

The abundances determined from the FCS fluxes and the CHIANTI atomic data seem 
to be in agreement with the standard solar value (An c /Ao = 0.15; Grevesse & Sauval 1998), 
even though the results from Figure 2 represent the highest possible value for the Ne/O 
abundance ratio for these active region spectra. These results are not in agreement with the 
Chandra values observed in stars, and therefore, do not support the hypothesis put forth by 
Antia & Basu (2005) - that an enhanced abundance ratio (A Ne /A = 0.52) would eliminate 
the discord between the solar models and the helioseismology observations. In light of these 
results, we have also reconsider several observations of the Ne IX and O VIII spectral lines 
from the pioneering days of X-ray astronomy. These full-Sun instruments had the advantage 
of observing the Sun as a star, as Chandra does, and were therefore not weighted by the 
smaller FCS field of view. 

(1) Blake et al. (1965) describe results from a U.S. Naval Research Laboratory Aerobee 
rocket flight on 1963 July 25. The payload included Bragg crystal spectrometers, pinhole 
cameras, and slit scanners. The spectral lines detected were attributed to a single active 
region near disk center. These lines included Ne IX A13.44, O VIII A16.01, and O VIII A18.97 
with fluxes of 2.1e4, 5.0e4, and 3.6e5 photons cm -2 s -1 , respectively. Using the CHIANTI 
atomic data and a temperature of Log T = 6.3, as we did for the FCS data, we find a 
maximum value for the neon-to-oxygen average abundance ratio of A^ e j Aq = 0.18. 

(2) Evans & Pounds (1968) observed the full Sun with two slitless Bragg spectrometers 
flown on a Skylark rocket on 1966 May 5. There were two overlapping sets of spectra 
from two visible active regions with intensities in ergs cm~ 2 s _1 of Ne IX A13.44 = 18., 7.8; 
Ne IX A13.70 = 7.4, 3.9; O VIII A15.18 = 12., 7.; O VIII A16.01 = 20., 10.; and O VIII A18.97 
= 130., 54. Averaging all the possible combinations gives A Ne /A = 0.25. 

(3) In a series of papers, Walker, Rugge & Weiss (1974) describe results from U.S. Air 
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Force satellites OV1-10 from 1967 January 4 and OV1-17 from 1969 March 20. The OV1- 
10 Bragg crystal (flat KAP) spectrometer and a proportional counter spectrometer which 
viewed the Sun directly were mounted on a compact solar pointer. The instrument detected 
intensities in ergs cm" 2 s" 1 for Ne IX A13.44 + A13.55 = 28.8; and O VIII A18.97 = 231. 
This gives A Ne /A Q = 0.20. 



3. Discussion 

The FCS data analyzed here and the rocket observations of the full X-ray Sun are 
consistent with the accepted solar abundance ratio of A^ e j Aq = 0.15. It is disappointing that 
these data do not confirm the prediction of Antia & Basu (2005) of A^ e /Ao = 0.52, but it is 
not necessarily surprising. We have no doubt that the harmony between the solar models and 
the helioseismology observations will eventually be restored. What is somewhat surprising, 
however, is that this abundance result is not consistent with the value of A Ne /A Q = 0.41 
obtained by Drake & Testa (2005) for their sample of 21 nearby stars observed with Chandra. 

There are, however, a few published results where the solar Ne/O abundance seems 
to be statistically greater than the accepted value. Murphy et al. (1991) analyzed data 
from the Solar Maximum Mission Gamma-Ray Spectrometer for the X5 limb flare of 1981 
April 27. They used a Bessel function for the accelerated particle spectrum and found 
A Ne /A > 0.35. Ramaty et al. (1996) used a power law rather than a Bessel function and 
revised the abundance to A^ e j Ao = 0.25. Schmelz (1993) studied FCS observations of the 
impulsive M5 flare on 1980 November 5 and found A Ne /A Q = 0.32. 

Since the Ne/O abundance ratio is the same for all layers of the solar atmosphere, these 
enhanced neon detections were explained via a mechanism that would preferentially dredge 
up neon from the photosphere. Shemi (1991) suggested that preflare soft X-ray radiation 
could penetrate deep into the solar atmosphere and create nonthermal ionization ratios at 
the base of the chromosphere. Because the photoionization cross section of neon is higher 
than that of oxygen, ionized neon could be selected along with the thermally ionized low-FIP 
elements for preferential transfer to higher levels of the solar atmosphere by the ion-neutral 
differentiating mechanism operating in the chromosphere. Therefore, the layer where the 
energetic particles and the ambient plasma interact during the flare would contain an over- 
abundance of ionized neon with respect to other high-FIP elements like oxygen. The model 
requires an extended period of energetic soft X-ray emission to build up the supply of ionized 
neon. This condition was met for both neon-enhanced flares: the 1981 April 27 flare was a 
long duration event which lasted several hours in soft X-rays and the 1980 November 5 flare 
had a long and intense preflare phase. 
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Since the Chandra spectra analyzed by Drake k Testa (2005) are most likely dominated 
by photons from very intense flares, it is possible that the enhanced neon abundances they 
measure result from a similar mechanism. Unfortunately, however, it does not seem likely 
that these higher neon abundances reflect the overall solar composition. If this result is true, 
than the proposed enhanced-neon hypothesis will not resolve the current controversy, and 
we will have to look elsewhere in order to restore the disrupted harmony between theory and 
observation. 

We would like to thank E. DeLuca, S. Basu, J. Drake, P. Testa, N. Brickhouse, and 
M. Coffee for inspiring discussions on abundances. JTS would like to thank the Harvard- 
Smithsonian Center for Astrophysics for their kind hospitality during her 2004-05 sabbatical. 
Solar physics research at the University of Memphis is supported by NSF ATM-0402729 and 
NASA NNG05GE68G. 
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Table 1. 



Ion 



A (A) Log T 



Transition 



Ne IX 
Ne IX 
Nc IX 
O VIII 
O VIII 
O VIII 
Fc XVIII 
Fe XVII 



13.45 
13.55 
13.70 
15.18 
16.00 
18.97 
14.20 
16.78 



6.6 
6.6 
6.6 
6.5 
6.5 
6.5 
6.8 
6.7 



Is 2 1 So - ls2s !Pi 
Is 2 1 S - ls2p 3 P! 
Is 2 1 S - ls2s 3 Si 
Is 2 S 1/2 -4p 2 P 3/2 ,l/2 
Is 2 S 1/2 - 3p 2 P 3/2 ,l/2 
Is 2 S 1/2 -2p 2 P 3 /2,l/2 



2s 2 2p 5 2 P 3/2 - 2p 4 (ld) 3d 2 D 5/2j3/2 
2p 6 1 S - 2p 5 3s 3 Pi 
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: Neon IX 13.55 A 

_= Oxygen VIII 18.97 A 

3.0x10 ' - 



r 2.0x10 13 
1.0X10 -13 




6.2 6.4 6-5 6.8 7.0 

Log Temp 



Fig. 1. — G(T) functions from CHIANTI for pairs of neon and oxygen lines (left) and ratios 
of these pairs (right), (a) Ne IX A13.45-to-0 VIII A16.01; (b) Ne IX A13.55-to-0 VIII A18.97; 
and (c) Ne IX A13.77-to-0 VIII A15.18. Each of the ratios has a local maximum at Log T 
= 6.3, the temperature used in Figure 2 to determine the maximum Ne/O abundance ratio 
consistent with the FCS active region spectra. 
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Fig. 2. — A/v e / 'Ao with F^ e /Fo from the FCS observations for the 20 spectra listed in Table 
2, and G (6.3)/GW C (6.3) from the CHIANTI database: (a) Ne IX A13.45-to-0 VIII A16.01; 
(b) Ne IX A13.55-to-0 VIII A18.97; and (c) Ne IX A13.77-to-0 VIII A15.18. The solid line 
is the weighted mean, dotted lines are ±la uncertainties, dot-dash line is A Ne /A = 0.41 
(Drake & Testa 2005), and the dash line is A Ne /A = 0.52 (Antia & Basu 2005). The 
weighted mean and its uncertainty are printed in the upper left corner of each plot. 
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Table 2. 



Date Time (UT) AR T(MK) NeA13.44 NeA13.55 NcA13.70 OA15.18 OA16.01 OA18. 



1 


1986 


May 20 


00 


:19- 


■01 


:12 


4729 


3.09 || 


4300.±320. 


1100.±300. 


3300 


±300. 


940.±550. 


4600.±660. 


35100.±3 


2 


1986 


May 21 


14 


:02- 


•14 


:55 


4731 


3.5r| 


3300.±300. 


1360.±260. 


2100. 


±270. 


1200.±740. 


5600.±750. 


49700.±4 


3 


1986 


May 23 


02 


:11- 


■03 


:07 


4731 


3-35' 


1850.±240. 


69O.±240. 


1760 


±240. 


1100. ±680. 


4200.±650. 


47600.±4 


4 


1986 


May 24 


04 


:54- 


■05 


:50 


4731 


3-47il 


3100.±300. 


300.±230. 


2400. 


±280. 


1100. ±720. 


4900.±650. 


46000.±5 


5 


1987 


Apr 11 


22 


:26- 


■23 


:21 


4787 


3.59;^ 


5300.±400. 


1200.±370. 


4400. 


±380. 


1600. ±800. 


55OO.±870. 


52400.±5 


6 


1987 


Apr 13 


01 


:08- 


■02 


:03 


4787 


3.80}| 


4300.±360. 


1000.±330. 


3700. 


±350. 


1700. ±750. 


4800.±680. 


53200.±4 


7 


1987 


Apr 14 


14 


:50- 


■15 


:45 


4787 


3.76 : °J 


4700.±340. 


700. ±300. 


3000. 


±370. 


1700.±800. 


5200.±670. 


61200.±4 


8 


1987 


Apr 15 


23 


:22- 


■24 


:17 


4787 


3.76 : }3 


7000.±690. 


1900.±600. 


5200. 


±600. 


2000.±1300. 


11300.±1900. 


64000. ±9 


9 


1987 


Apr 19 


14 


:41- 


■15 


:10 


4787 


3.85 : °| 


11000.±900. 


3200.±800. 


9400 


±800. 


3800. ±1500. 


12000.±2200. 


138000.± 


10 


1987 


Nov 27 


16 


:20- 


■17. 


:17 


4891 




3100.±330. 


68O.±260. 


2800 


±300. 


1100.±780. 


6100. ±1000. 


59000.±4 


11 


1987 


Nov 29 


20 


:07- 


■21 


:06 


4891 


3.47;il 


2400.±560. 


200. ±73. 


2700 


±600. 


2400.±1000. 


3100.±800. 


32300.±7 


12 


1987 


Dec 07 


14 


:48- 


■16 


:01 


4901 


3.63;1J 


3500.±300. 


720. ±230. 


2200. 


±290. 


1000. ±500. 


3400.±750. 


35600.±4 


13 


1987 


Dec 08 


09 


:49- 


■10 


:51 


4901 


4.03;i| 


5100.±400. 


1100.±320. 


3400 


±350. 


1300.±620. 


4200.±890. 


40700. ±4 


14 


1987 


Dec 09 


03 


:06- 


■04 


:08 


4901 


3.94-0§ 


3500.±370. 


740.±300. 


2200. 


±320. 


1600.±640. 


5100.±770. 


48500.±4 


15 


1987 


Dec 11 


02 


:13- 


■03 


:13 


4901 


3.8o;°i 


4900.±350. 


820.±300. 


3000 


±310. 


1700. ±780. 


6400.±970. 


60600.±4 


16 


1987 


Dec 11 


10 


:04- 


■11 


;04 


4901 




3700.±360. 


900.±300. 


2500 


±320. 


740.±690. 


5300.±1000. 


41500. ±4 


17 


1987 


Dec 13 


09 


:10- 


■10 


:09 


SELimb 


3-31;i| 


6400.±440. 


1400.±410. 


6100. 


±430. 


3100.±1030. 


10000.±850. 


86200.±5 


18 


1987 


Dec 15 


03 


:3d- 


■04 


:31 


4906 


3.72;}f 


5900.±420. 


1100.±380. 


4200 


±380. 


2300.±800. 


6800.±790. 


71000.±5 


19 


1987 


Dec 16 


07 


:49- 


■08 


:46 


4906 


3.80-J! 


4300.±360. 


620.±200. 


3600. 


±420. 


1300. ±730. 


5400.±740. 


55800.±4 


20 


1987 


Dec 18 


03 


:46- 


■04 


:42 


4906 


2.95$ 


1700.±270. 


56.±19. 


1000. 


±230. 


870.±520. 


3100.±600. 


30600.±4 



